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Abstract 

A method used recently to obtain a discrete formalism for classical fields with nonlocal 
actions preserving chiral symmetry and uniqueness of fermion fields yields a discrete 
version of Huygens' principle with free discrete propagators that recover their continuum 
forms in certain limit. 
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1 Introduction 



Lattice theory has been used since many years ago as one of the non-perturbative ap- 
proaches to study physical effects that could occur in QCD or QED |T|. Nevertheless, the 
discrete formalism yielded by the Lattice theory does not preserve important properties 
such as chirality and uniqueness which are present in the standard field theory. Nielsen 
and Ninomiya @ have shown that for a local and translationally invariant hermitian dis- 
cretization of the fields it is not possible to have simultaneously chiral symmetry and 
uniqueness. Many attempts have been made to circumvent this restriction. The use of 
nonlocal operators can be found among others (see for example in spite that the 

locality of the Dirac operator is in general a desired property. In this context, a new ap- 
proach to the problem of discretization of fields is presented in J7J. This approach is based 
on a non-equispaced lattice, performed through the zeros of Hermite polynomials, where 
both discrete derivatives and Fourier transform have support || |{|. The technique yields a 
projection of the quantum algebras on a finite linear space yielding matrix representations 
for the partial derivatives that produce discrete hermitian actions with nonlocal kinetic 
terms. The fermion doubling and the chiral symmetry breaking are absent in such a for- 
mulation. In this work we obtain a discrete version of Huygens' principle for free spinor 
fields in 1 + 1 D. We show that whenever the number of the nodes in the non-equispaced 
lattice tends to infinity, the propagator we found approaches to their continuum form. 
We also show that the discrete propagator is indeed a discrete Green function by taking 
the inverse of the Dirac operator. 

This work is organized as follows. In section |2| we review the essentials of our discretiza- 
tion method and introduce a discrete version of the sign and step functions. In section 
| we obtain the discrete Huygens' principle on the non-equispaced lattice. Section (| is 
devoted to study the discrete version of the free Dirac propagator in 1 + 1 D. Finally in 
section | we give our main conclusions. 



2 Discrete technique 
2.1 Review of the method 

In this subsection we only present the main results of our discretization scheme; proofs 
and further applications can be found in [0, ^ [10], [11], ^ . 

Let us consider the non-equispaced four dimensional lattice constructed with the set of 
nodal points Xj (fi denotes the Lorentz index and j = 1,2,..., N^, with being the 
number of nodes along the direction /i) performed by the zeros of the Hermite polynomial 
Hn^O- We denote N = N NxN 2 N 3 as the number of lattice points. From this set of 
points we construct the four N x N matrices: 

D = D <g) 1^ ® 1jv 2 ® ljvi, D 3 = ljv ®.D3®ljv2®ljvi, (1) 
D 2 = l No ® 1jv 3 ®D 2 (g) l Nl , D x = l No ® 1jv 3 <g> l N , 2 ® Di, 
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where 1^ is the identity matrix of dimension and is the skew-symmetric matrix 

0, i=j, 



jk 



a u ? 

rf r rr*' 



i^3- 



Let us define the matrix S = So (g> S3 (g> S2 <8> Si where (S^)^ = 5jk exp[—(xj) 2 /2]H' Nii (xj). 
Then, SD^S -1 is a projection of the partial derivative <9 M in the subspace of functions 
U generated by products of the form u n (£) = exp(— £ 2 /2)if n (£), n = 0, 1, . . . , — 1, 
with £ = x^ and /i = 0, 1,2,3. This means that such matrices are exact representations 
of the partial derivatives for functions in U. Therefore, whenever iV M — > 00 we will get 
convergent approximations to the partial derivatives of a function ip(x) spanned by the 
basis {wn(0}o°- Thus, the discrete version of a dynamical differential variable operating 
on such functions is essentially the matrix obtained under the replacement <9 M — > SD M S _1 . 
The error arising from this procedure can be estimated in special cases JTJ] and can be 



related to the complement of ip(x) with respect to U M. Thus, if ip(x) G U and \I/ denotes 
the Af x 1 vector of components 

* q = il)(x q )=$(x),3? h ,x z u x Q m ), (2) 

ordered according to q = j + (k — l)Ni + (I — l)NiN 2 + (m — l)NiN 2 N 3 where first j runs 
over j = 1, . . . , Ni, then k over k = 1, . . . , N 2 , then / over 1 = 1,..., N 3 , and finally we 
take m = 1, . . . , N . In this form we have that the vector ^ constructed with the values 
of ip )ft = d^ip at the site x q = (x^, x\,x\, xj), is given by 

= SD^S- 1 ^. (3) 

The similarity transformation given by S changes \I/ into itself except for an alternating 
change of sign along each direction when — ► 00 J/], therefore we may use D M instead 
SD^S" 1 as a discrete representation of <9 M . 

Let g(x fl ) be a given function and g and g' the vectors of components g(xj) and dg(xj)/dx^, 
respectively. Let us denote by G = diag(g) the diagonal matrix whose nonzero elements 
are the components of g. Then, we have that 

(S^D^S^ )G = G(S ll D fl S~ 1 ) + G' + R^, 

where G' = diag(g') and R is the residual matrix which projects an arbitrary vector on 
the orthogonal subspace generated by the complement basis {w n (£)}^ . This equation, 
applied to the vector constructed with the values of the function h(x^) at the nodes, is 
the finite representation of the familiar formula d(gh)/dx^ = gdh/dx^ 1 + hdg/dx^. Since 
the complement basis becomes empty as goes to infinite, we have that 

D„G = GD^ + G', -> 00, (4) 

where G' can be substituted by the diagonal matrix diag(_D M g). 

In summary, D M becomes a discrete representation of yielding convergent results for 
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a wide class of functions. However, this operator is not a local one from the field theory 
point of view [Q, ffl, i.e., there exist nonzero elements of D M for which 

KD/Jgg/I > exp(-\\x q - XfWoo) (5) 

since l/\z\ > exp(— \z\). Here || • is the max vector norm. In Sec. || we will see a 

consequence of this property. 

Now, let us define the symmetric function 

N-l 

F{t,ri) = EWVKOwfa). (6) 

1=0 

where ^(0 = (2 N - 1 ~ l {N - l)!/(JVn)) 1/2 ^,(0/^-i(0- Since 

iD F = FP°, -iDj-F = FP J , (7) 

the commuting matrices (|l|) can be diagonalized simultaneously by the unitary and sym- 
metric matrix 

F = f£®F 3 ®F 2 ®F 1 , (8) 

where (F^jj. = F(xj,p%) is a discrete Fourier transform for one variable |9]]. Here, pj 
is also a zero of Hjy (£) and it represents an eigenvalue of the discretized momentum. 
The elements of the four-dimensional discrete Fourier transform (§) satisfy the important 
asymptotic formula 

lim F q , q = C N e-^'V, (9) 

iV — >oo 

where p q ■ x q / = Pm x m' ~ Pf x f' ~ v\ x \> ~ p] x ]' anc ^ Cn is the product of the constants 
C Nfi = 2 Ar "~ 3 / 2 (r[(iV / ,+l)/2]) 2 /jV At ! except for an alternating change of sign g [7|. By using 
an asymptotic expression for the Hermite zeros we find that Cn becomes proportional 
to the standard measure of a Riemann integral in each variable (the difference between 
two consecutive lattice points): 

CW„ = 4=Ax" = 4=A^, (10) 



where Ap M = Ax^ 1 = 7r/y2A^. The equations (0) and ( [TOD are the backbone of the 
quadrature formula for the Fourier transform ||. 
Additional and useful properties of F M , are 

= F li U = UF lt , (11) 

where C/ is the AT M x matrix whose entries are given by = 5j t N„-k+i- Denoting by 
[■}k the /cth column of a matrix, the first equality means that 

ftj*- (-l) N » +1 [Fi]k, (12) 

under p£ — >■ — p%, p = 0, 1, 2, 3. 
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2.2 Discrete distributions 



In section f| a representation of the step function 0(t — t') on the lattice will be needed to 
construct the discrete propagators. In order to do that, let us give a discrete form of the 
sign function e(t — t'). Is it well known that 



e(t - t') = / = / e ltE — e~ lt E dE. 

nJ-oo E H J-oo E 



According to (|9]) and (|1(J) the integral with the integrand e ltE becomes the matrix \^2~kF 
applied to the remaining part of the integrand. Thus, the discrete form of e(t — t') is the 
skew-symmetric matrix 

where P° is the diagonal matrix whose nonzero elements are E^. Of course, 5 has to 
be a nonsingular matrix, restricting us to consider Nq even. The fact that the entries 
of H are real is guaranteed by (0), and ©-(0) yields a finite asymptotic value for any 
fixed element^] of H. Since D Fq = —iF^P , the discrete derivative of H is two times the 
identity matrix divided by the measure, i.e., 

A,H = 2(^)l Wo . (13) 

The measure goes to zero as N — > oo, therefore, this result is the discrete version of 
de{t — t')/dt = 25 (t — £'). ^From ( |13"1) we see that -iFoiP )-^ is the inverse of Do and 
that D E - ED = 0. 

On the other hand, the step function accepts the Fourier representation 



i r°° e^-^dE _ i r°° e^'^dE e /■*> e i{t -^ E dE 
t[ ' ~ ' ' ~ 2^ J-oo E-ie ~ ~2^ J-oo E + 2^ J-oo ~E 2 ' ^ ' 



This equation is just the simple relation 6{t) = [e(t) + l]/2. Since the first positive zero 
of Hn (E) (the first positive value of Ek, i.e, Enq_ .) goes to zero as N — > oo, we may 

2 + 

give a discrete representation of 6{t — t') through (O). The identification of e with En^ 

2 +i 

gives the matrix 

= i(S + E), (15) 
as the representation of 6{t — if) on the lattice, where we have defined 



AE 



By a fixed element of a, N x N matrix A, we mean the entry A N / 2 +j,N/2+k, with j and k fixed. Here 
N is even and j, k — ±1, ±2, . . . , ±N/2. 
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A fixed element of this symmetric matrix becomes asymptotically equal to 1 divided by 
the measure, i.e., 

= ^ W - oo. (16) 

Note that the discretization of (|T4]): 

and (|15D are equal up to the first order in En$_ . The application of D to E yields 
—Enq E. Since a fixed element of 5 is bounded as N — ► oo and Eno — > 0, we get the 

2 ' 2 ' 

expected property 

A)e = 5, (17) 

where 5 is the identity matrix divided by the measure, i.e., 5 = 1n xn / 'AE. In order to 
normalize the right-hand side of (|17|) is convenient to our purpose define 

6 = -*F (P° - ^^W^r^o = \{ 2 + S) (18) 

as the discrete form of 0(t — £'), where 5 = (A.E)E and E = (A.E')E. In this way we have 
that -Do© = Ia^oxVo- Another important relation between these matrices follows from (^). 
If we choose the vector g as the jth column of O, i.e. [Q]j, Eq. (fD becomes 

A)diag([eL-) = diag([e] J )L> + [W^li, iV -> oo, (19) 



3 Huygens' Principle 

The above formalism is applied in Ref. J7| to find discrete spinor fields and some of their 
properties. In particular, the discrete version of the free Dirac equation is found to be 

iY <g> = mU N ^, (20) 

where 7 M are Dirac matrices, I^n is the identity matrix of dimension 4N and \I> is 
the discretized field whose four spinorial components \l/ a have spatial-temporal indexes 
(^ a )q = ^a(x q ) ordered according to @. We remark that the Dirac operator D = ry^tgiD^ 
is not a local one. To show this, let us take the matrix norm induced by the max vector 
norm 
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\\A NxN \\ 00 = N max \a jk \. 

l<j,k<N J 

Thus, according to (|5[), we have that there exist nonzero elements of D for which 

||D(xg, x 9 lloo > max |(D M ) gg /| > exp(-\\x q - vlU)- (21) 

0</^<3 
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Accordingly, the no-go theorem due to Nielsen and Ninomiya does not apply in this case. 
Since the the condition for chiral symmetry is satisfied trivially by D H, this operator 
yields chiral fermions with no doublers. 

The square of a mass eigenvalue of (|2DD, say m r , is given by the discrete form of the 
energy-momentum relation 

m 2 r = {plf ~ ip)? ~ (Pi) 2 ~ (P?) 2 , (22) 

where pf is a nonzero component of the diagonal matrix representing the discretized 
momentum p^. Such values are zeros of Hermite polynomials and therefore, m r is degen- 
erate. The indexes in (|22"D are ordered according to 

r — j + (k- 1)N! + 1)N!N 2 + (to - l)JViiV 2 iV 3 + (i - 1)N^N 2 N Z N Q 

where the slowest index is i = 1,2,3,4. The degeneracy of m r does not depend on 
and the discrete mass-shell consists only in few points even when — > oo. However, 
since the Hermite zeros become dense on the axes, mf, approaches to any real number 
as Np — > oo and its degeneracy is given by the number of points sufficiently close to the 
mass-shell defined by m r . This behavior is illustrated in Fig. 1. 
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Figure 1: The degeneracy of m r is approximated by the number of points pf of (22) 
located between the hyperbolas defined by \m — m r \ = e/2 corresponding to the 
dashed lines in this figure for the case m r = 7.3. The points satisfying this condition 
are only displayed and the distance between the hyperbolas has been exaggerated. 
Since we can always find points pf close enough to any point laying on the lower or 
the upper hyperbola, both surfaces are completely covered when Nu — ► oo. 



Thus, the sum over the discrete mass-shell can be substituted asymptotically by a sum 
over all the values p\, k = 1,2,3, i = 1,2, ... , Nk, and p° given by 

(P°J 2 ~ {p)f - (Pi) 2 - (Pf) 2 = ™ 2 , (23) 
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where m is now a parameter whose numerical value is close to m r . This approximation 
will be used later. 

On the other hand, a general discrete solution of (|2Cf ) can be written as a linear com- 
bination in the degenerate subspace corresponding to the eigenvalue m r . By using (|T^), 
this combination can be written as a sum over positive values of the energy and other 
sum over the negative values. Since the spinor space is not affected by the discretization 
procedure we can obtain the positive and negative energy bispinors following the standard 
procedure. Thus, a plane-wave solution of (|20|) corresponding to the mass eigenvalue m r 
is the vector 



[Ft] 



(24) 



where al and b\ L Z 



1,2, are complex numbers and 
u!, = C mr [l, 0,pf/(p° m + m r ),pt/(p° m + m r )f, 



u 



C mr [0, l,p fc /(p^ + m r ), -pf/(p° m + m r )f 



3 

Cmr\Pi/(Pm + ™r)>Pj/(Pm + ™r), 1, 0] T , 



Here, C mr = (wir/Pm)[(Pm + m r )/2m r ] 1//2 (the w spinors are normalized to one), pf k = 
p x - ± ip^ and j, k, I, m are the indexes of the components of p 2 q satisfying (p2|) for a given 
eigenvalue m r . Eq. (p2]) is the condition that defines the sums given in (|24T): positive 
energies on the discrete mass-shell correspond to the first sum and negative energies to 
the second one. Since C mr depends on the inverse of p^ (eigenvalue of P°) the number 
of temporal nodes must be an even integer to have pJJ, / 0. 

To obtain a relation between the temporal components of \l/ mr we need to write down 
the explicit dependence on their spatial and temporal indexes maintaining the spinor 
structure; thus, the order of the tensor product of ([24|) should be taken according to this. 
We begin by taking only 1 + 1 variables to illustrate a standard procedure that can be 
generalized to more variables. The representation of the Dirac matrices employed here is 



a x and 7 1 = ioy. Taking into account flT2"|), Eq. (^) can be written as 



m r )q 



^m r (xj,t k ) = E + a(pj',E k ')u{p j r,E k t)F*(t k ,E k >)F(x j ,p jl 



+b(-p jf , -E k/ )v(-p f , -E k r)F*(t k , -E k r)F(xj, -py 



(25) 



where the sign appearing in (|T2| ) has been included in the complex constant b q > 
b(—pj/, —E k >) and the relations between indexes are 



where j,j f = 1, 



q = j + (k 
N% and k, k' = 



1)N, 



Q r — f + (k r — l)iVi, 
A^ . Here, the spinors are given by 



m r / sfE + p 
y/E~+p 



-P, 



-E) = — 
; 2E 



VE^p 



-m r / y/E - p 



(26) 
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If we multiply (p5|) by F*(xj,Pj„) and sum over Xj we obtain 

^2F*(x j ,p j „) ty mr (xj,t k ) = J2 + a(pj',Ek>)u{pf,E k t)F*(t k ,E k/ )S. 



P\,= m l 



PjlhPjl 



+b(-p f , -E k ,)v(-p f , -E k/ )F*(t k , -E k ,)5 Pjll _ P] , 
where we have made use of the orthogonality relations 

J2F*(xj,pj„)F (xj ,pf) = S jf/jl = 5 p .„ tP .,, 

Xj 

ZF*(x j ,p jfl )F(x jl -p r ) = (Jjv/.jvi-j'+i = 5 Pj „,- P l . 



(27) 



The right-hand side of (|27|) is different from zero if pj„ is on the discrete mass-shell. By 
using the orthogonality relations for the spinors and projectors we obtain that 



^m r (xj,t k ) = ^2G(xj,t k ;x j ',t k >)j 0li if mr (x j/ ,t k/ ), 



[28) 



where 



G(xj,t k ; Xj',t k >) 



v\„=m? r 2-Efc// 



' ^ q n + rn r ) iy + j ' ,V J'\ F\t k , E k „)F(xj,pj„) 



+0/ 



qll 



F* (ty , E k „ 
F*(xj>, -p jff 



F*{t k/ ,-E, 



kii) 



F*(t k , -E k// )F(xj, -pj„) 



(29) 



Eq. (|28|) seems to be a discrete version of Huygens' principle but a simple calculation 
shows that it is just a resemblance. 

It is true that such an equation yields an identity for ^ mr (xj, t k ) as given by (pop, when 
k — ► (interchange the sums of (p8|) and sum first over Xf), but the function given by 
does not satisfy G(xj,t k i]Xj>,t k ') = 7°^^' and can not be considered as a discrete 
Green function. 

However, according to the argument given at the beginning of this section, we can sub- 
stitute asymptotically the sum over the discrete mass-shell given in (^9|) by the sum over 
all the values of the discrete spatial momentum satisfying a energy-momentum relation 
according to (p3|), i.e, 

£ + - E (so) 



Ni, k = 1,...,N . Here, m is a parameter close to m r 



as — ► oo for j = 1, 
Therefore, in this limit (p9|) can be substituted by 



2K 



k'l 



i" + mrJ F*(t k ,,E k/l) 



-F*(t k ,E k „)F(xj,p j „) 



mr] F*{xf, P^ F *( tk , -E k „)F(xj, - Pj „) 



F*(t k i, —E k „) 
9 



(31) 



It is not difficult to show that 

G(xj,t k >;xj>,t k >) = 7%',j', 
so that we can interpret ( |31~D as the discrete Green function and 

^ mr (Xj,h) =J2G(Xj,h;Xji,tkih ^rn r (Xj>,tk>), ty < t k , (32) 

x j' 

as the discrete version of the Huygens' principle. The condition t k i < t k has been added 
to define the retarded discrete propagator. 



4 Propagators on the lattice 

Let us denote by K(xj,t k ;xf,t k i) the 1 + 1 retarded propagator. Thus 

K(xj,t k ;Xj',t k i) = G{Xj,t k ;Xj',t k '), t k i < t k . 



(33) 



In this section we will obtain some properties of this function. First, we will show the 
equivalence between ([33]) and the discrete propagator given in |§. By using the asymptotic 
formula fl9T) we have that 



lim J-(o,aQ = lim F*(o,- Pr) 

iv M ^oo F*(0, E k >) n^oo F*{0, -E v ) 



Therefore, 
K(x,t;0,0) 



E 

y/p 2 +m 2 =E 



1 

2E 



'+m)F*{t,E)F(x,p) + (p / -m)F*{t,-E)F{x,-p) , (34) 



where t > and the indexes have been suppressed to clarify the notation. Writing 
^(£> v) = v) + ^(^j if) an d using the parity of these functions with respect to their 
arguments, we get the components 



im 



K n (x,t;0,0) = K 22 (x,t;0,0) = £ - —S(t,E)C(x,p) 



■\J p 2 +m 2 =E 



K 12 {x,t; 0,0) = V C{t,E)C(x,p) - |s(t, E)S(x,p) 



\/ p 2 +m 2 =E 



K 21 (x, t; 0, 0) = Yl E ) C ( X ,P) + f^(*, E)S(x,p) 



\/ p 2 +m 2 =E 



E 



which are the same that those given in [g], obtained here under a more general scheme. 
As it is shown in fl, this discrete Green function converges to the correct continuum 
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propagator and can be rewritten as the weighted sum over lattice paths of the checkerboard 
model of Feynman. 

As usual, the condition t k i < t k in (|33|) can be included as the product of the step and 
Green functions 

K(xj,t k ; Xf,t k ') = Q(t k , t k t) G(x j} t k ; Xf,t k >), (35) 

where is defined by (|18|). The Feynman propagator can be defined along the same lines. 
The Hadamard product of two matrices A and B [defined as (A o B)j k = Aj k Bj k ] can be 
used to write (|35"D in matrix form. Let us define the (2iV) x (2N) matrix 

= l 2x2 ® 6 ® O, 

where N = N Nx and O is the identity matrix for the Hadamard product (Oj k = 1, 
j, k = 1, 2, . . . , Ni). Thus, taking the order for the tensor product used to obtain ([H]), we 
have that 



Now, if F denotes the matrix of components (F) qq > = F(xj,pj/)/ F(t k , E k r), the matrix 
block G(xj,t k ;xf,t k r) given by ( pi]) takes the compact form 



G 



qq' 



E 

v /p 2 j „+m 2 =E k „ 



(F)^(- 



- hm ^(F%,q + (F% q ^ 



E, 



kll 



qll 



m r 



E, 



kll 



Iq/lq 



(36) 



and (pBD becomes the product (&) qq i(G) qq i. This means that the retarded propagator is 
the (2N) x (2N) matrix 

K = 0oG. 



By using the relation 

(A® B) o (C ® D) = (Ao C) ® (B o D), 
it is not difficult to see that the components of K can be written as 



(K) 



qq' 



E 

^P 2 J „+m 2 =E k/l 



e o fJ) 



F x 



iq„ + 



9 o F ) fi 



qll 



w" E k „ 
— m r 



i q"<i 



+ 



qq" E, 



kll 



iq"q 



Applying (27^ ® D M — m r l 2 Ar) to (|37| ) and using (0), (|T9| ) and (0) we get 



2ArJgs(K) S(? 



gg'- L 2x2 



(37) 



(3? 



(39) 



as — >• 00. Thus, the retarded propagator defined by (|37D becomes asymptotically equal 
to —i times the inverse of the Dirac matrix: 



K 



(40) 
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Similarly, the Feynman propagator S which propagate the positive frequencies toward 
positive times and the negative ones backward in time, is defined by 



Iqq' 



- 2i 



v /p2 // +m 2 =E k „ 



e o f$) ® f 1 



qll 



+ m r 



qll 



11" E k „ 
— m r 



i qi'q 



qqll 



E, 



kn 



I qllq 



(41) 



Both S and K are related as in the continuum case. By using ([18|) and taking into 
account that 5 and S are skew-symmetric and symmetric matrices respectively, we have 
that 6 T = -6 + S. On the other hand, flT6|) means that £ = O as iV — > oo, therefore 



/<29 



+ 



E 



(F f ; 



99" 



m r 



v /p i» +m2 



few 



The second term of the right-hand side of this equation is a solution of (fy M ® D^ t — 
^rl2Af)^ — 0. Since K satisfies (^9|), S becomes the inverse matrix of the Dirac operator. 
We end this section by discussing the form that such a matrix should have for a finite N^. 
Strictly speaking, there is no inverse for (ij^ ® D M — m r l 2 Ar) since m r is an eigenvalue 
of ® D^. However, according to the approximation procedure given at the beginning 
of section [| we can substitute m r by a value m ^ m r laying close enough to m r to yield 
a nonsingular Dirac operator. As a matter of fact, this value plays just the role of an 
e-prescription in the propagator to avoid singularities on the real axis. 
Proceeding in this way and taking the representation for the 7-matrices used before, we 
obtain 

-1 _ f mA AD_ 



S = (i^ ® D M - ml 



(42) 



where the N x N matrices A and D± are defined as follows 



fD D 



m 2 l N )' 



D± = i(D ±Di). 



(43) 



Since D and Di commutes, 



Nj 



so that A is the inverse of the discrete representation of the Klein-Gordon operator 
—6q + df — m 2 . It is not difficult to see that A, D M and D± also commute. Thus, ([42]) 
takes the more compact form 



S = {ij" ® D M + ml N ){l 2x2 A) 



(44) 



which is the discrete form of the well-known relationship between fermionic and scalar 
propagators. Note that this result does not depend on any asymptotic limit; we only need 
to take N% 7^ -/V and m 7^ m r to have nonsingular finite matrices. 
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5 Final remarks 



We end this work by pointing out the main differences between the discrete technique for 
a finite and its asymptotic formulation given in section [| 

First of all, we would like to remark the similarity and convergence of the results of this 
method to the ones of standard field theory. As it was shown in sections |^ and f|, the 
convergence of the discrete Fourier transform it is not enough to obtain a propagator as 
a sum of plane-wave solutions of the equations; it is needed to include other terms laying 
in the neighborhood of the mass-shell. This is also a requirement in the standard field 
theory where the propagators are assumed off-shield and Eq. is no longer satisfied. 
On the other hand, in a physical process such as dispersion of particles, the energy and 
momentum conservation are consequence of translational invariance which is accounted 
by Huygens' principle. In our discretization scheme the translational invariance is lost; 
however later is recovered when the number of the lattice points goes to infinity, since the 
Hermite zeros become dense in the real axis. 

Finally, we have shown at the end of the previous section how to get a finite matrix 
form of the propagators without reference to any asymptotic limit. This matrices can be 
computed easily in numerical calculations of this discrete field theory. 
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